Reconstructing Fourier's law from disorder in quantum wires 
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The theory of open quantum systems is used to study the local temperature and heat currents 
in metallic nanowires connected to leads at different temperatures. We show that for ballistic wires 
the local temperature is almost uniform along the wire and Fourier's law is invalid. By gradually 
increasing disorder, a uniform temperature gradient ensues inside the wire and the thermal current 
linearly relates to this local temperature gradient, in agreement with Fourier's law. Finally, we 
demonstrate that while disorder is responsible for the onset of Fourier's law, the non-equilibrium 
energy distribution function is determined solely by the heat baths. 

PACS numbers: 72.15. Jf,73.63.Rt,65.80.+n 



The search for a microscopic derivation of Fourier's law 
[1], or even a microscopic demonstration of it, is still a 
major theoretical challenge [2, 3]. In recent years, the 
problem seems more relevant than ever due to both the 
continuous miniaturization of electronic devices and the 
need for alternative energy sources. Both these trends re- 
quire better understanding of fundamental processes of 
energy transport in nanoscale systems. Following, there 
have been many theoretical attempts to derive Fourier's 
law in various systems, both classical and quantum [4- 
10]. In the quantum regime, attention has been focused 
mainly on small spin-chains [4, 7, 11] or quantum har- 
monic oscillators [12], where it was demonstrated that 
Fourier's law holds for chaotic systems, but vanishes 
when transport turns ballistic. This is in agreement with 
the behavior of classical systems, where it has long been 
postulated that chaos leads to Fourier's law [4]. 

However, little is known on the energy transport in 
nanoscale electronic quantum systems. One of the rea- 
sons for this is that in order to demonstrate the validity 
- or violation - of Fourier's law one needs to (i) define 
a local temperature out of equilibrium and show that it 
develops a uniform gradient, (ii) evaluate the local heat 
current, and (iii) show proportionality between these two 
quantities. The first task is especially difficult since tem- 
perature is a global equilibrium property, and it is not 
clear if a " local temperature" can be defined at all when 
the system is out of equilibrium [13, 14]. For this reason, 
recent studies of the origin of Fourier's law either use a 
phenomenological definition of local temperature (as an 
expectation value of a local energy operator [4, 15]) or as- 
sume that a temperature gradient is already present [16]. 
An alternative route is to study the energy diffusion in 
closed systems (i.e. without thermal baths) [17]. 

An additional reason that renders calculation of energy 
transport in electronic systems a formidable task is the 
fact that the size of the Hilbert space scales exponentially 
with the number of electrons, making numerical calcula- 
tions very demanding. This is why previous numerical 
calculations on heat transport in quantum systems usu- 
ally refer to very small systems, typically of the order of 
ten spins (see, however, Ref. [18]). 



Here we report a calculation of energy transport in 
electronic quantum wires that overcomes both the above 
issues. It is based on solving the quantum master equa- 
tion for non-interacting electrons in the presence of dis- 
sipative baths (held at different temperatures) in the 
Markov approximation. We use a recently suggested 
method [19, 20] to map the many-electron problem to a 
single-particle system, which allows calculations for sys- 
tems an order of magnitude larger than in previous stud- 
ies, and enables us to define a local temperature oper- 
ationally, that is one directly accessible experimentally, 
even out of equilibrium. 

Model. - The system consists of a linear chain bonded to 
small leads, which are connected to thermal baths held at 
different temperatures (upper panel of Fig. 1). [19, 20] 
A similar set-up was used for spin-chains in, e.g., Ref. 
[21]. The Hamiltonian of the system is given by H = 

H L +H R + Hd + H c , where Hl,r,4 = 12ieL,R,d e i c \ c i ~ 
tJ2{ij)£~L Rd ( c 1 c j + h.c.^j arc the tight-binding Hamil- 
tonians of the left lead (L), right lead (R) and wire 
(d, of length Ld), respectively (t is the hopping inte- 
gral, which serves as the energy scale hereafter). H c = 
(gLC L Cd,o + gnc^Cd.Ld + h.c.) describes the coupling be- 
tween the left (right) lead to the wire, are creation 
operators for an electron at the point of contact between 
the left (right) lead and the wire, and Cd,o {cd,L d ) de- 
stroys an electron at the left-most (right-most) sites of 
the wire. The on-site energies ti are randomly drawn 
from a uniform distribution U[— W/2, W/2], with W be- 
ing the disorder strength. The lattice constant is taken 
to be a = 1, and we consider here spinless electrons. 

The quantity of interest, from which the required in- 
formation (such as local temperature, density and heat 
current) may be extracted, is the single-particle den- 
sity matrix, defined by p — ~^2 kk i Pkk'\k){k'\ , where 
p kk ' = Tv{c\c w pmb) , pmb is the full many-body 
density-matrix, and |fc) are the single-particle states. In 
Ref. [19, 20] it was shown that p obeys a master-equation 
of the Lindblad form [22] (setting h = 1) 



-t[H,p} + C L [p} + C R [p} = C[p], 



(1) 
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where £l(r) [p] are super-operators acting on the den- 
sity matrix, describing the left (right) thermal baths, 
held at temperature T L ( R \, and in contact with the left 
(right)-most side of the leads (solid lines in the up- 
per panel of Fig. 1). The super-operators are defined 
in the Lindblad form [22, 23] via V-operators, C[p] = 

Efe.fc' (-H V kk> V kk>,p} + Vkk>pv£ k ,y with {•,•} being 
the anti-commutator. The T^-operators are generalized to 
account for the different baths, and are given by [20, 21] 

L,K) 
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where fn' R \ek) 



1/ (cxp ^ ^bTlr ) + ^) are tne Fermi distributions of the 
left and right leads, with \i the chemical potential. The 
coefficients 7^,'^ = J2 ri es L , R V'fcfa) 7o ipk'( T i) (where 
tpk (r) are the single-particle wave functions) describe the 
overlap between the single-particle states and \k') 
over the region of contact Sl(r) between the left (right) 
baths and the corresponding junction leads, shown by 
the solid lines in the upper panel of Fig. 1. 70 describes 
the strength of electron-phonon (bath) interaction. The 
form above can be derived from first principles by tracing 
out the bath degrees of freedom, with the latter formed 
by a dense spectrum of boson excitations (e.g., phonons), 
which interact locally with electrons at the edges of the 
system. 

Local temperature. - The non-equilibrium steady state 
of the system is given by the solution of the equation 
£[/?]= 0. In order to calculate the local temperature, 
we attach a third thermal bath (described by an addi- 
tional term £ t i p (r)[p] in the master equation) which is 
connected to a given position r of the wire and serves as 
an external probe (upper panel of Fig. 1). The resulting 
equation for the steady state with inclusion of this probe 
now reads (L + £ t i p (r))[p] = . The temperature T t ; p 
of the additional operator is changed, and the resulting 
local density, n(r) = ^2 k Pkk\fpk(r)\ 2 [24], is compared to 
that obtained without the additional operator. We then 
define the local temperature as that T t ; p for which there 
is minimal change in the local density. [20] Physically, 
this corresponds to placing a local temperature probe in 
close proximity to the wire [25]. When the probe has the 
same local temperature of the wire, there is no heat flow 
between them and thus the local properties of the wire 
arc unchanged. 

Figure 1 shows the local temperature profile along a 
chain of length Ld = 150. The numerical parameters are: 
lead dimensions L x — L y = 3, temperature of the left and 
right leads is TL = 0.1 and Tr = 1 respectively, and elec- 
tron number tie — 56 (corresponding to one third filling). 
The disorder strength is W = (solid line), W = 0.1 
(small dashing) and W = 0.5 (large dashing) and aver- 
aged over 500 realizations of disorder (all calculations in 
this report were performed over a wide range of param- 
eters, yielding similar conclusions). One can see several 
features from Figure 1. The most prominent feature is 



the fact that for clean (W — 0) and weakly disordered 
(W = 0.1) wires the temperature hardly changes along 
the wire. This is clear evidence that Fourier's law is in 
fact violated under these conditions. Instead, a uniform 
temperature gradient at the center of the wire is recov- 
ered for large disorder (W = 0.5). In addition, from 
Fig. 1 one can see an asymmetry between low and high 
temperatures - the local temperature of the wire never 
reaches that of the colder bath (with Tl = 0.1), but 
does reach the hotter bath temperature (with Tr = 1). 
This asymmetry can be attributed to the structure of the 
Fermi function, which gives different weights to high and 
low temperatures. Also notice the thermal length at the 
edges of the wire, where the local temperature is roughly 
constant. 

In order to determine a relation between the appear- 
ance of a uniform temperature gradient and the onset 
of chaos [4], the inset of Fig. 1 shows the distribu- 
tion function P(s) for the level-spacing s of neighboring 
single-particle energy levels, for W = 0.1 (thin line) and 
W = 0.5 (thick line). The dashed line is the Wigner- 
Dyson distribution, Pwd(s) = sexp(— as 2 ) (where a is 
some constant), which is conjectured to correspond to the 
onset of quantum chaos [26] . As seen, the distribution for 
weak disorder shows features reminiscent of the ordered 
system and is very different from the Wigner-Dyson dis- 
tribution. Accordingly Fourier's law is invalid. On the 
other hand, for strong disorder we see both agreement 
with the Wigner-Dyson distribution and a uniform tem- 
perature gradient at the center of the wire, suggesting 
that indeed these have the same physical origin. 

Heat current. - The next step in validating Fourier's 
law comes from calculating the local heat current j(r), 
and evaluating the thermal conductivity n such that 
j(r) = — «(r)VT(r) . It is clear that for a uniform tem- 
perature profile (as that of a clean wire), a finite heat 
current results in a divergent k. On the other hand, for 
the system presented here one expects that the heat cur- 
rent is always finite, since energy is always injected and 
extracted from the system. It is thus natural to define 



a global thermal conductivity, K 



Ldj 



where j is 



t r -t l ■ 

the heat current averaged over the whole wire. Follow- 
ing this definition, it is evident that K = R (where R is 
the average over the wire of the local thermal conductiv- 
ity, defined via Fourier's law) only when Fourier's law is 
valid. This definition thus gives us a mean to point to 
the onset of Fourier's law. 

To make this argument substantial one must calculate 
the local heat current. Since there is some ambiguity 
in defining the local heat current, we follow Ref. [16, 27] 
and define the heat current operator via a time-derivative 
of the local energy operator Hi, i.e. = — ji 

(in units of t 2 , taking K=l). The local energy opera- 
tor in the wire is defined as Hi — ti\i)(i\ — ^t(\i)(i + 
1| + — 1| + H.c). The time-derivative of a general 
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FIG. 1: (color online) Local temperature as a function of po- 
sition along the wire. The local temperature is calculated for 
three different values of disorder, W = (solid), 0.1 (small 
dashing), and 0.5 (large dashing). For clean wires (W = 0) 
the temperature is uniform, and a uniform temperature gra- 
dient builds up as the disorder increases. Upper panel: ge- 
ometry of the model system. The solid lines at the edges 
correspond to the contact area of the thermal baths. Inset: 
energy level spacing distribution for W = 0.1 (thin line) and 
W = 0.5 (thick line). For weak disorder the distribution 
shows a structure reminiscent of the clean system. For strong 
disorder the distribution resembles the Wigner-Dyson distri- 
bution (dashed line), marking the onset of chaos. 



operator O in the Lindblad formalism is given by [23] 

^ = i[H,d] + E kk ,(-H v L>Vkk>>d} + v kk ,dvZ k ,). 

From this relation it is straightforward to find ji and 
to calculate its expectation value, ji. Fig. 2(a) shows 
ji as a function of position along the wire (same pa- 
rameters as in Fig. 1), for different values of disorder 
W = 0, 0.1, 0.5. Note that in the central region of the 
wire the heat current is uniform, and increases at a dis- 
tance of the order of the thermal length of the leads (solid 
lines in the W = 0.5 curve are guides to the eye), due to 
the proximity of the heat baths. The inset of Fig. 2(b) 
shows the averaged heat current (over the whole wire) 
as a function of temperature difference AT — Tr — Tl, 
for different values of disorder. From the linear regime 
one can extract the global K, shown by empty circles in 
Fig. 2(b) as a function of disorder. On the other hand, 
one can extract the averaged local k by using the local 
heat current (Fig. 2(a)) and the local temperature (from 
Fig. 1), calculated and averaged over sites close to the 
center of the wire. The values of the local k are shown 
as filled circles in Fig. 2(b), and exhibit a divergence for 
clean wires. Only at large disorder, where Fourier's law 
is valid, do the two definitions of thermal conductivity 
coincide. In fact, we can take this as criterion for the on- 
set of Fourier's law, which microscopically corresponds 




1 t 

20 40 60 80 100 120 140 
position 




0.0 0.1 0.2 0.3 0.4 0.5 



W[t] 

FIG. 2: (color online) Local heat current and thermal con- 
ductivity, (a) Local heat current j as a function of position 
for different disorder strengths, W = 0,0.10.5. The current 
is largest close to the leads, where the contact with the heat 
baths is located. The solid line on W=0.5 is a guide to the 
eye. (b) Local (solid circles) and global (empty circles) ther- 
mal conductivity as a function of disorder. While the global 
thermal conductivity hardly changes with disorder, the local 
thermal conductivity shows strong disorder dependence, and 
diverges in the clean limit due to the vanishing temperature 
gradient. The two tend to the same value only with increasing 
disorder, when Fourier's law becomes valid. Inset: averaged 
heat current as a function of temperature difference AT for 
different disorder strengths. From the linear regime of these 
curves the global thermal conductivity was evaluated. 



to the value of disorder W of order t/2, namely when 
the strength of the disorder is comparable to half of the 
kinetic energy parameter. 

In order to further determine the role of disorder, we 
study the local and global electron energy distribution 
function, both of which can be measured experimentally 
[28] . In terms of the density matrix, the global distribu- 
tion function is simply /(-Efc) — p^k, and the local dis- 
tribution function is given by f\ oc {E k ,r) = \ip k (r)\ 2 p kk . 
Fig. 3 shows the distribution function for the clean 
wire (W = 0, squares) and disordered wire (W = 0.5, 
circles, same parameters as in Fig. 1) as a function 
of energy. As seen, there is hardly a difference be- 
tween the two functions. The solid line is the curve 
f(E) = \{h(E) + f R (E)) , where f L{R) (E) is the Fermi 
distribution of the left (right) lead, with the correspond- 
ing temperature. The excellent agreement between the 
numerical curves and the averaged density suggests that 
one cannot consider the wire as a system with an effective 
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FIG. 3: (color online) Global electron distribution functions 
with and without disorder. The distribution function is eval- 
uated for a clean wire (squares) and disordered wire, W = 0.5 
(circles). The thick line corresponds to an average distribu- 
tion f(E) = ±(/l(£) + ME)) , where f HR) (E) is the Fermi 
distribution of the left (right) lead, with the corresponding 
temperature. Inset: the local electron distribution function 
at the edge of the wire (thick line) and the center of the wire 
(thin line) for a clean system (upper panel) and a disordered 
system (lower panel). The resemblance between the two cases 
demonstrates that the distribution function is determined not 
by the local electronic structure, as dictated by the Hamilto- 
nian, but rather by the boundary conditions, i.e. the thermal 
baths. 



temperature T e g = |(Tl+Tr.) ,but rather it is the distri- 
bution function itself that is averaged, and that disorder 
does not affect the distribution function. 

In the inset of Fig. 3 we plot the local distribution 
function for the clean and disordered wires at the left 
edge of the wire (thick line) and at the center of the wire 
(thin line). Again, the curves for clean and disordered 
wires are very similar. We note that in our geometry, in 
contrast to mesoscopic wires [28] , one cannot attribute to 
the local distribution a simple position-dependence, even 
in the presence of disorder. 

These findings show that the energy distribution func- 
tion is only determined by the boundary conditions, 
namely those provided by the bath operators, and not by 
the local structure (i.e. the Hamiltonian) . On the other 
hand, nature of the heat transport is determined solely 
by the microscopic character of the Hamiltonian states. 
This implies that Fourier's law cannot be validated from 
measuring the local distribution function. 

One possible way to verify Fourier's law is to study 
the local temperature. The local temperature of elec- 
tronic systems can be measured experimentally [25]. Rel- 
evant experimental systems for which our predictions 
may be tested can be, e.g., carbon nanotubes [29], quan- 
tum point contacts [30], atomic-size metallic wires [31] or 
silicon nanowires [32, 33], where measurements of ther- 
mal conductance and thermo-power have already been 
demonstrated. A different route is to measure the length- 
dependence of the thermal conductivity, as was recently 



applied to nanotubes [34] 

Finally we note that the above theory did not include 
the effects of electron interactions. While for some quasi- 
one-dimensional systems electron interactions may not be 
very important, it is nevertheless of interest to find out 
whether electron interactions alone suffice for the valid- 
ity of Fourier's law, even in clean wires. However, the 
above theory is limited to non- interacting systems, and 
to resolve this issue one must employ a more elaborated 
method that can encompass both interactions and envi- 
ronments [35] . Such studies are currently underway. 

This work is supported by DOE under grant DE-FG02- 
05ER46204. 
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